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ABSTRACT 

A general three-field variational principle is obtained for the motion of an acoustic fluid enclosed in 
a rigid or flexible container by the method of canonical decomposition applied to a modified form 
of the wave equation in the displacement potential. The general principle is specialized to a mixed 
two-field principle that contains the fluid displacement potential and pressure as independent fields. 
Semidiscrete finite-element equations of motion based on this principle are displayed. 

1. INTRODUCTION 

An elastic container (the structure) is totally or partly filled with a compressible liquid or gas (the 
fluid ) . The fluid structure system is initially in static equilibrium in a steady body force field such 
as gravity or centrifugal forces. We consider small departures from equilibrium that result in forced 
or free vibratory motions. To analyze these motions the fluid is treated as a linear acoustic fluid, 
i.e., compressible but irrotational and inviscid. The purpose of the present work is 

1. To derive variational equations of motion based on a mixed variational principle for the fluid 
subsystem. 

2. To obtain semidiscrete equation of motions following spatial discretization of the coupled prob- 
lem by the finite element method. 

The derivation of the mixed variational principle for the fluid is based on the method of canonical 
equations advocated by Oden and Reddy [l] for mechanical applications. The most general dynam- 
ical principle derived in this paper contains three primary variables: pressure-momentum vector, 
dilatation-velocity vector, and displacement potential. 

The general principle is specialized to a two-field functional of Reissner type that has pressure 
and displacement potential as primary variables, as well as a free coefficient a that parametrizes the 
application of the divergence theorem. The coupled variational equations are then discretized by the 
finite element method. Linkage with the elastic container is made to establish coupled semidiscrete 
equations of motion for a rigid container and for the fluid-structure system. By appropriate selection 
of the coefficient a a number of discrete formulations of the elastoacoustic problem results. The 
computational implications of these selections are briefly discussed. 
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Table 1 Notation for Fluid States 


Quantities 

Domain 

Original 

Reference 

Current 

Transient 

Displacements 

V 

0 

d° 

d* 

d = d*-d° 

Velocities 

V 

0 

d° 

d‘ 

d=d'-d° 

Boundary displacements* 

S 

0 

d° n 

d*n 

<U = d‘-d« 

Displacement potential 

V 

0 

ij> 0 



Pressures 

V 

0 

P° 

P* 

p = ,'- p o 

Body forces 

V 

0 

b = V/? 

> 

II 


Density 

V 

p 

P 

P 


* Positive along outward normal 


2. GOVERNING EQUATIONS 

The three-dimensional volume domain occupied by the fluid is denoted by V. This volume is 
assumed to be simply connected. The fluid boundary S consists generally of two portions 

S : Sd U S p . (1) 

S d is the interface with the container at which the normal displacement d n is prescribed (or found as 
part of the coupled fluid-structure problem) whereas S p is the “free surface” at which the pressure 
p is prescribed (or found as part of the “fluid slosh” problem). If the fluid is fully enclosed by the 
container, as is necessarily the case for a gas, then S p is missing and S = S d . The domain is referred 
to a Cartesian coordinate system (xj, ® 2 , z$) grouped in vector x. 

The fluid is under a body force field b which is assumed to be the gradient of a time independent 
potential /?(x), i.e. b = V/?. All displacements are taken to be infinitesimal and thus the fluid 
density p is invariant. 

We consider three states or configurations: original, from which displacements, pressures and 
forces are measured, current , where the fluid is in dynamic equilibrium at time t, and reference, which 
is obtained in the static equilibrium limit of slow motions. Transient motions are the difference 
between current and reference states. It should be noted that in many situations the original 
configuration is not physically attainable. Table 1 summarizes the notation used in relation to these 
states. 

2.1 Field Equations 

The governing equations of the acoustic fluid are the momentum, state and continuity equations. 
They are stated below for the current configuration, and specialized to the reference configuration 
later. The momentum (balance) equation expresses Newton’s second law for a fluid particle: 

pd* = — Vp* + b = -Vp* + V/?. (2) 

The continuity equation may be combined with the linearized equation of state to produce the 
constitutive equation that expresses the small compressibility of a liquid: 
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p* = —KVd* = — pc 2 Vd‘, (3) 

where K is the bulk modulus and c = y/K/ p the fluid sound speed. If the fluid is incompressible, 
K,c — + oo. This relation is also applicable to nonlinear elastic fluids such as gases undergoing 
small excursions from the reference state, if the constitutive equation is linearized there so that 
K = po(dp/dp)o. 

The boundary conditions are 

d*=d* on S d , P*-P* on S p , (4) 

where d„ is either prescribed or comes from the solution of an auxiliary problem as in fluid-structure 
interaction, and p may be either prescribed or a function of d n and b, as in the surface-wave (“slosh”) 
problem. 

2.2 Integral Abbreviations 

In the sequel the following abbreviations for the volume and surface integrals are used: 

{f)v d = [ fdV , [g] s d = f gdS, (y] s<£ d = f gdS , etc. (5) 

Jv Js Js d 

That is, domain-subscripted parentheses (square brackets) are used to abbreviate volume (surface) 
integrals. Abbreviations for function innerproducts are illustrated by 

(f.g)v =' / fgdv, U,g)vxt d = [“ [ fgdVdt, (/,j] S jX . =' f' / fgdSdt, etc. 
Jv Jto Jv Jto Js d 

( 6 ) 

3. THE DISPLACEMENT POTENTIAL 


S.l The Reference State 

Taking the curl of both sides of (2) yields 


curld* = 0. (7) 

The general integral of this equation for a simply connected domain is 

d* = V0* + a + bt, (8) 

where rp* = t/>‘(x,f) is the displacement potential, a = a(x) and b = b(x) are time-independent 
vector functions, and t denotes the time. If accelerationless motions (for example, rigid body 

motions) are precluded by the boundary conditions, a and b vanish. Replacing d = V0* into the 
momentum equation (2) we get 

Vp* = -pVrp* + V/?, (9) 

which spatially integrated gives 

p* = -pfr + /? + C{t), (10) 
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where the scalar C{t) is not spatially dependent. Next integrate the constitutive equation (3) over 
V and apply the divergence theorem to Vd: 


( P*)v + (pc 2 Vd*)v = (p 1 )^ + [pc 2 d„] s — 0. (11) 

Inserting p 4 from (10) into the above equation furnishes a condition on C(t), which gives 

c(t) = -^-K\s + ^‘)y - i(/J)v = -^-K\s + pT‘ -7, (12) 

where v = (l) v is the fluid volume and / = (/)v/v denotes the volume average of a function / 
defined over V. Substituting C(t) into (10) we get 

P* = -p(^ 4 - 0«) + (/?-/?)- ~[d^\ s . (13) 

In the static limit of very slow motions, the inertia terms may be neglected and we recover the 
reference solution 

p° = (/?-?)-vW ]*• ( 14 ) 

For an incompressible fluid \d n ]s = 0 but c — ► oo; thus it would be incorrect to conclude that 
p° = /? — /?. A counterexample to this effect is provided in [2]. 

S.2 Transient Motions 

Subtracting the constitutive relations at the current and reference states we get 

p = — pc 2 V 2 V> = pc 2 s , (15) 

where s = — V 2 V> is called, following Lamb [3], the condensation. Subtracting (14) from (13) yields 

P = -P{4> ~ VO - ~[d n \s- (16) 

On equating (15) and (16) we get modified forms of the wave equation that account for mean 
boundary surface motions: 


s » = t-jL + or = (17) 

C XJ 

The second form follows from —vs = [d„.]s, which is a consequence of the divergence theorem. For 
an incompressible fluid, c — ► 00 and [d n ]s = 0, and from the first of (17) we recover the Laplace 
equation V 2 0 = 0. 
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S.S Adjusting the Displacement Potential 

If the transient displacement potential is modified by a function of time: 

0 = 0 + P(t), (18) 

where rp is the potential of (8)-(17), we may chose P(t) so that c 2 0 = V 2 0 = —s for any t. Then 
we obtain the classical wave equation 

= or (—5 - c 2 V 2 ) = 0. (19) 

In the sequel it is assumed that this adjustment has been made. If so, C(t) vanishes and (16) 
reduces to 

P = ~pi>- ( 20 ) 

4. MIXED VARIATIONAL PRINCIPLES 
4.1 Canonical Decomposition 

In this section we derive multifield variational principles for the fluid domain following the canonical 
decomposition method advocated by Oden and Reddy [1], This method is applicable to self-adjoint 
boundary value problems (BVP) of the form 


Au — f in D (21) 

where u is the unknown function, / the data, A a symmetric linear operator, and D the domain of 
existence of the solution. For time-dependent problems D is the tensor product of the time domain 
(typically 0 to t) and the volume V. To apply this method, the operator A is factored as 

Au — T*ETu = /, (22) 

where T and E are linear operators in V and T* is the adjoint of T. This is called a canonical 
decomposition . This decomposition may be represented as the operator composition sequence 

Tu = e, Ee = a, T*a = /, (23) 

where e and o denote intermediate field variables in D. The three equations (23) are called the kine- 
matic, constitutive and balance equations, respectively, in mechanical applications. The canonical 
representation of boundary conditions on the surface S = 5 U U S„ is 

Bsu s = g on S u , Bgcr s = h on S (24) 

where Bs and B* s are surface operators, g and h denote boundary data, and ug = ls u a ^d 
<*s = T go are extensions of u and o to the boundary S. The extension operators 75 and 6g often 
involve normal derivatives. 


5 



4.2 The Wave Equation 

The classical wave equation (19) is not a good basis for the canonical decomposition (22). Its 
principal drawback is that the pressure field does not appear naturally as an intermediate variable 
in (23). A more convenient form for our purposes is obtained by taking the Laplacian of both sides 
of (19), and multiplying through by the density p : 

pV 2 (0 — c 2 V 2 0) = 0, whence A = pV 2 t / = 0. (25) 


A suitable canonical decomposition is A = TET t , where 


T = 



E = p 


I 0 

c 5 


0 - 2 



(26) 


in which * = y/—l. Boldface symbols axe used for T and E because these are 4x1 and 4x4 
matrices, respectively. The operator product sequence (23) becomes 


(27) 

The intermediate fields e and a are 4x1 column vectors. These vectors are partitioned into their 
temporal and spatial derivative subvectors for convenience in subsequent manipulations. Note that 
the transient pressure p appears naturally as the spatial component of a. The temporal components 
of e and a are the complex velocity iv and complex specific momentum tm, respectively. 

The boundary portions S u and S 9 of (24) are relabeled S d and Spy respectively, to match the 
notation (1). Boundary and initial conditions may be stated in generic form as 


Brj)(x,t) = g(x,f) on S dt BV(x,t) = h(x,t) on S p , d(x,t 0 ) = d 0 (x), d(x,ti) = d^x). 

(28) 

Here B and B* are time-independent 4x1 and 1x4 vectors, respectively, related to the canonical 
B s and B| operators of (24) by B = B sis and B* = B^Ts, where 7 s (a scalar) and Ts (a 4 x 4 
matrix) are boundary extension operators for rp and <r, respectively. Comparison with (4) and the 
use of Green’s function reveals that 

Bf = -B* s = [0 0 0 1 ], g T = [0 0 0 <f n ], ls = ^~, r 5 = I, h = -p. (29) 


4-S Three Field Principle 

The most general variational principle for the canonical decomposition (26) allows the three fields: 
0, e, and <r, to be varied independently. The principle may be stated as 6L(tp,e,a) = 0, where the 
functional L is [2] 

L(u,e,<r) = Ly+Ls = §(Ee,e)v*t + (<r,TV>- v)v*t - (/, rp)vxt + (o^B^-g^xt - (h, rp s )s r xt, 

(30) 
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where Ly and Ls collect volume and surface terms, respectively. On inserting (27-29) into (30) we 
get 


L v = A(Ee,e) V xt + (*,T0 -e) Vxt = £ j ^ [|p(~v T v + c 2 s 2 ) - m T (V0 - v) - p(V 2 0 + s)] dV d 


L s = {*s, B0 - g)$ d xt - ( h,xff S )s,xt = 



dt. 


whereas (/, 0)vxt vanishes. 
4-4 Two Field Principles 


(31) 


A two field principle of Reissner type can be derived from the functional L by enforcing the inverse 
constitutive equations e = E -1 <r a priori. The resulting principle, which allows 0 and o to be 
varied simultaneously, is 5J2(0, <r) =0, where 


R(0,<r) = R v +R S = - 5 (E~ 1 <r,<r)vxt + (<r,T0)v'xt-(/.0)vxt + (*s,B0-g)5 (J xt-(>i,0s)5,xt. 

(32) 

where Rs = Ls and 


Rv{*l>,<r) = -|(E _1 <r,(r)vxt + (<T,T0)vxt= f \ (^-m T m - - m T V0 - pV 2 0) dV dt. 

Jto Jv 2 p 2 pc 

(33) 

The specific momentum disappears as an independent field if we enforce m = pV0 a priori , 
whereupon the functional R becomes a function of 0 and p only and the volume term contracts to 


Rv(*,p) = 1' J v (~ heWf?* ~kfz~ P v V) iV <“• 

To check R — Ry{4),p) + Rs we form its first variation: 

6R = - (pV 2 ,? + V 2 p,S^j v - (-ijp + V 2 <l>,Sp) v + [pV 2 ,j + Vp,6<l>] 

-[p-P.J^l -[f*-4»*p], -[ 

I dnls r xt lan Js^x* L anJt 0 


(34) 


Sxt 


(35) 


Setting 6R = 0 provides the field equations and boundary conditions. 

4.S Parametrization 

A one parameter family of variational principles can be obtained by transforming all or part of 
the last term in (34), viz. pV 2 0, by the divergence theorem (Green’s first formula for the Laplace 
operator) 



xJjdV + 


[ V0VpdF= f p^ds = 
Jv Js dn 



dxb 

>Tn dS + 



(36) 
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Let 0 < a < 1 be the portion of that term to be transformed. Insert pV 2 ip = apV 2 ip + (l — a)pV 2 ip 
in (35) and apply (36) to apV 2 ip to get 


R a v = J* 1 y v (ip(VV>) T V^-i^+aVpVV'-( l-«) P V 2 il>)dV-aJ s Pj^dS-a |£, pdS 

(37) 

Finally, replace the Laplacian V 2 ip left over in (37) by c -2 (0) to arrive at the parametrized two-field 
functional 


dt. 


Ratty, P) = RaV +Rs = £ f v {~ W V ^) TV ^ ” 2^2 + aV P V ^ “ (l “ <*)c ~ 2 pV>) ^ 




(38) 


dt. 


The highest spatial derivative index for both primary variables ip and p is 1, except if a = 0, in 
which case it is only 0 for p . The two interesting limit cases are of course a = 0 and a = 1, for 
which 


^-s:[L (-*#(v* r v^-i§-.-v)^+ j s ,{?±-i)is+ j S '%pds 

*i(<M = /' [/,(- -^ + vpvv-) rfv - ^ pd„<is| dt. 


C It, 
39) 
(40) 


5. FINITE ELEMENT DISCRETIZATION 


5.1 Discretization of R a 

In the following we derive semidiscrete finite-element equations of motion based on the R a functional 
(38). The volume V is subdivided into fluid finite elements. Over each fluid element the state is 
represented by the primary variables ip and p, which are defined as functions of position in the usual 
shape-function interpolation procedure. The finite element interpolation in V may be expressed as 

ip{x,t) = Nj(x) 9(t), p(x,t) = Nj(x) p(t), (41) 

where 9 and p are computational column vectors that contain nodal values of ip and p, respec- 
tively, and N^, and N p are corresponding arrays of dimensionless shape functions. The specified 
displacement over Sd is interpolated by 

d»(x, t ) = n T d(x, t) = n r Nj (x)d, (42) 

where n is the external-normal unit vector on Sd, N d contains the displacement shape functions of 
the enclosing container, and d contains nodal displacement values. For the moment these displace- 
ments will be assumed to be prescribed, hence the superposed tilde. Inserting (4l)-(42) into the 
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functional (38) yields the semidiscrete quadratic form 


J2 a (¥,p) = -|p^ T H*-^-p T Gp + a® T Fp + (l-a)[® T Wp-® T Dp + « f %] -p T T r d, (43) 


where 

H = [ VNjVN^, dV = H t , F = [ VN^VN 4,dV, G = f c~ 2 NJN p dV = G r , 

Jv Jv Jv 

I) = J <T 2 NjNp<iV, W =J (V„N*) T Np«i5, T t = J Nj n T Nj<iS, f * = J g f V » N * <iS - 

(44) 

The integration with respect to time is dropped as it has no effect on the variation process described 
below. 

5.2 Continuity Requirements 

The interelement continuity requirements of the shape functions of xp and p depend on the index of 
the highest spatial derivatives that appears in R a . If a ^ 0, this index is 1 for both xp and p and 
consequently C° continuity is required. It is then natural to take the same shape functions for both 
variables: 

N* s N p (45) 

with both vectors $ and p of equal dimension and evaluated at the same nodes. Then' some of the 
matrices in (44) coalesce as 

H = F, G = D. (46) 

The case a = 0 is exceptional in that no spatial derivatives of p appear. One can then chose C~ l 

(discontinuous) pressure shape functions; for example, constant over each fluid element. If this is 
done, obviously 

N* # N p (47) 

because xp must be C° continuous. Furthermore the dimensions of p and will not be generally 
the same. 

5.S The Rigid-Container Equations of Motion 

Since R a contains time derivatives of of order up to 2 in 9, the appropriate Euler-Lagrange varia- 
tional equation is 


6R a 


dR a 

d9 


d dRa 
dt d9 


a 2 dRg\ 
dt 2 89 / 


69 + 


dRa 

dp 


6p=0, 


which applied to (43) yields 


[pH» + aFp - (1 - a) Dp + (1 - a) Wp + (1 - a)fy] 69 = 0, 
[ - p -1 Gp + aF T 9 - (1 - a)D r ¥ + (1 - a) W r ¥ - T r d] 6p = 0. 


(48) 


(49) 
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These equations can be represented in partitioned matrix form as 


pH 

[—(1 — a)D T 


VKtHA -AolCM 


-(1 - o)f* 1 
T r d J ’ 


(SO) 


where J — (1 — a) W + aF. 

5.4 The Flexible-Container Equations of Motion 

If the fluid is enclosed in a flexible container, the boundary displacements d are no longer prescribed 
on Sd but must be incorporated in the problem by including them on the left hand side of the 
equations of motion. We shall only consider here the case in which the container is modelled as a 
linear undamped structure for which the standard mass /stiffness semidiscrete equation of motion is 

Md + Kd = t d + Tp, (51) 


where M is the mass matrix, K the tangent stiffness matrix at the reference state, Tp is the pressure 
force on the structure, and f d is the externally applied force on the structure. Note that K in general 
must account for prestress effects through the geometric stiffness. Combining (50) and (51) we get 
the coupled system 


M 

0 

0 


f d ] 


' K 0 

-T 


r d 1 


f U } 

0 

pH 

-(1 - a)D 

< 

1 * 1 

► + 

0 0 

J 

i 

' a# 1 

► = 4 

-(1 - a)U 

0 

-(1 - a)D r 

0 


Ip J 


-T t J t 

, 

0 

1 

1 


Ip J 


0 . 


(52) 


! If a = 0, 


MOO 
0 pH D 
0 D t 0 



(53) 


There is little than can be done beyond this point, as the shape functions for p and will be generally 
different. Although the pressure may be constant over each element, no static condensation of p is 
possible in the dynamic case. 


| 

2 

O 

O 

1 


d ] 

K 0 -T 


fd) 

\U 

O 

H 

0 

i 

* \ + 

OOF 

< 


0 > 

[ 0 00 


Ip J 

-t t f t -p~ x g 


Ip J 1 

l 0 J 


( 54 ) 


Note that all these systems are symmetric. 


5.5 Identical Shape Functions 

Further progress in the case a = 1 can be made if we assume, as discussed in §5.2, that the shape 
functions for p and t/> coincide. Taking then (46) into account, (54) simplifies to 


fM 0 O’ 


d] 

K 0 -T 


fd] 


{fa) 

, 

O O 

W 0 
0 0 


: * 

0 0 H 

— T r H —p~ l G 

< 

Or 

Lp J 

► = < 

° 

l® J 


( 55 ) 
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The second matrix equation gives pH\5 r + Hp = 0. Since H is nonnegative definite we must have 


P = 


-p4f. 


(56) 


This is the discrete analog of the continuous relation (20) for the dynamic overpressure. It is possible 
to rearrange (54) into several equivalent symmetric forms; for example, 


' M pT O' 

f d 1 


K 0 0 

pT T -pH G 


► + 

0 0 0 

0 G 0 

l P J 


0 0 p -1 G 




(57) 


plus two more symmetric forms not shown here. Note that (55) and (57) are not equivalent for the 
static limit, since the static matrix of the latter becomes singular. In fact, only (55) handles that 
condition in a numerically stable way. 

5.6 Static Condensation 

If (56) is used to eliminate the pressure vector from (55) we obtain 


[“ ?]{}}♦[-* !]{i}-{S)- 

Conversely, eliminating the displacement potential vector gives 

[# :h;h; ■» k;}-m 


Unlike the previous equations, both (58) and (59) are unsymmetric. Thus the elimination of a field 
variable, be it p or ip, causes symmetry to be lost. 

6. CONCLUDING REMARKS 

Displacement-potential formulations are of practical interest in fluid-structure transient-response 
and vibration problems since they form the basis for effective numerical computations. For some 
recent applications see [4-10] and references therein. The preceding treatment unifies a number 
of previous continuum-based and algebraic statements [7-12] of the coupled problem. It may be 
extended in the following directions: 

(1) Elastoacoustic free vibrations of the coupled system. The symmetry properties of the algebraic 
eigenproblem are studied in [8-13]. 

(2) Gravity waves and free-surface “slosh” modes. An additional field, the elevation r/, appears 
over S p ; but at the same time the fluid may be taken as incompressible. See e.g. [8]. 

(3) The inhomogeneous wave equation c 2 V 2 ^» — ip = /, / ^ 0, when the body force field b(x ,t) is 
time-dependent and V 2 b ^ 0. Additional forcing terms, including f p , appear in the equations 
of motion. These are of interest for rotating fluids. 

(4) Retaining the specific momentum m as independent field. 
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SUMMARY 

The equations of motion for an acoustic fluid enclosed in a moving or flexible container 
are studied. It is shown that the determination of the reference state must account for the 
surface-integrated effect of the wall motions. The governing equation of transient motions 
about this state in the displacement potential does not generally reduce to the classical 
wave equation unless special adjustments are made. The results are relevant to finite 
element formulations based on the displacement potential. 
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1 Problem Description 

The results presented in this Note were obtained in the course of a wider study by Felippa 
and Ohayon (1988) of variational methods for transient motions and vibrations of acoustic 
fluids held in flexible and/or moving containers. These partial results merit special atten- 
tion on two counts. First, they extend the classical dynamic equations of acoustic fluids 
to include wall motions as well as the static limit in a consistent manner. Second, they 
are relevant to finite element implementations that have not accounted for the correction 
terms described herein. Computations that are particularly affected by these corrections 
involve liquid masses subject to prescribed dynamic motions, such as tanks and reservoirs 
under seismic ground-motion excitations, and rocket fuel tanks under launch conditions. 

The general problem is as follows. A container (the structure) is totally or partly 
filled with a compressible, homogeneous liquid or gas (the fluid). Although the container 
is generally flexible, the rigid-but-moving container case is not excluded. The fluid is 
modelled as an acoustic medium (the linearly-compressible generalization of an ideal fluid) . 
We consider dynamic motions about a static reference state, which will be determined as 
part of the study. 

If the container is rigid and fixed, the reference state is well known: the static equi- 
librium solution in which the pressure is equal to the hydrostatic pressure, and the dis- 
placements may be taken as zero. The acoustic motions about this position are governed 
by the homogeneous wave equation in the displacement potential. But if the container 
walls move, we show that a correction term that depends on the mean boundary motion 
appears. The reference state is affected, and the resulting transient vibration problem is 
no longer given by the classical wave equation unless special adjustments are made. 

A boundary-integral term representing the mean container motion was introduced by 
Aganovic (1981) for the surface-wave problem of an incompressible fluid posed in terms 
of the velocity potential. Ohayon (1987) considered similar terms in the displacement 
potential formulation. The general forms presented in this Note for a compressible are 
believed to be new. 

2 The Acoustic Fluid 

The three-dimensional volume occupied by the fluid is denoted by V. This volume is 
assumed to be simply connected. The fluid boundary S consists generally of two portions 

S : S d U S p . (1) 
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Table 1 Notation for Fluid States 


Quantities 

Domain 

Original 

Reference 

Current 

Transient 

Displacements 

V 

0 

d° 

d* 

d = d*-d° 

Velocities 

V 


d° 

d* 

d = d*-d° 

Boundary displacements* 

S 


dn° 

d' n 

d n =d*~ d « 

Displacement potential 

V 



V 

Ip Ip* — Ip® 

Pressures 

V 

0 

P° 

P* 

w -0 
P-P -P 

Body forces 

V 


b = V/3 

b = V/3 


Density 

V 

9 

9 

9 


• Positive along outward normal 


Sd is the interface with the container at which the normal displacement d n is prescribed 
(or found as part of the coupled fluid-structure problem) whereas S p is the “free surface” 
at which the pressure p is prescribed (or found as part of the “slosh” problem) . If the fluid 
is fully enclosed by the container, as is necessarily the case for a gas, then S p is missing 
and 5 = Sj. The domain is referred to a Cartesian coordinate system (xi, X2, X3) grouped 
in vector x. 

The fluid is under a body force field b which is assumed to be the gradient of a time 
independent potential /?(x), i.e. b = V/3. All displacements are taken to be infinitesimal 
and thus the fluid density p is invariant. 

We consider three states or configurations: original , from which displacements, pres- 
sures and forces are measured, current , where the fluid is in dynamic equilibrium at time t, 
and reference , which is obtained in the static equilibrium limit of slow motions. Transient 
motions are the difference between current and reference states. It should be noted that in 
many situations the original configuration is not physically attainable. Table 1 summarizes 
the notation used in relation to these states. 

Field Equations. 

The governing equations of the acoustic fluid are the momentum, state and continuity 
equations. They are stated below for the current configuration, and specialized to the ref- 
erence configuration later. The momentum (balance) equation expresses Newton’s second 
law for a fluid particle: 

pd* = -Vp t + b = -Vp t + V/?. (2) 

The continuity equation may be combined with the linearized equation of state to produce 
the constitutive equation that expresses the small compressibility of a liquid: 
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p 4 = -JCVd* = -pc 2 Vd 4 , (3) 

where K is the bulk modulus and c = \jKfp the fluid sound speed. If the fluid is 
incompressible, K,c —*■ oo. This relation is also applicable to nonlinear elastic fluids such 
as gases undergoing small excursions from the reference state, if the constitutive equation 
is linearized there so that K = Po{dp/dp)o. 

The boundary conditions are 

d* = i 4 on S d , p 4 = p 4 on 5 P , (4) 

where d£ is either prescribed or comes from the solution of an auxiliary problem as in 
fluid-structure interaction, and p may be either prescribed or a function of d n and b, as in 
the surface-wave ( ‘‘slosh*) problem. 

3 The Reference State 

Taking the curl of both sides of (2) yields 

curl d * = 0. (5) 

The general integral of this equation for a simply connected domain is 

d 4 = Vifr* + a + bt, (6) 

where V ,< = V' t ( x >0 Is the displacement potential, a = a(x) and b = b(x) are time- 
independent vector functions, and t denotes the time. If accelerationless motions (for 
example, rigid body motions) are precluded by the boundary conditions, a and b vanish. 
Replacing d =Vrp t into the momentum equation (2) we get 

Vp 4 = -pVV>* + V/?, (7) 

which spatially integrated gives 

p* = -pt* +0 + C(t), (8) 

where the scalar C(t) is not spatially dependent. Next integrate the constitutive equation 
(3) over V and apply the divergence theorem to Vd: 

(p‘) v + d‘V = (p‘)v + [pSd'Js = o. (9) 

Inserting p 4 from (8) into the above equation furnishes a condition on C(t), giving 

c(«) = -^-KU + £ = ~K\s + p¥-p, (io) 
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Figure 1. Cylindrical fluid container in gravity field 


where v = (l) v is the fluid volume and / = [f)v/v denotes the volume average of a 
function / defined over V. Substituting C(t) into (9) we get 

j>* = (ii) 

In the static limit the inertia terms may be neglected and we recover the reference solution 

p° = 03 - ?) - ^(<1*. (12) 

For an incompressible fluid, [d„.]s = 0 but c — ► oo; thus it would be incorrect to conclude 
that p° = 0 — 0. To illustrate this point* consider a rigid cylindrical container of cross 
section area A, filled with liquid up to height H — H\ + H 2 . The origin of the Cartesian 
system (x x , x 2 , X3) is placed at H 2 below the free surface, with x x = x upwards and 
normal to that surface; see Figure 1. The body force is the gravity field b = (— py,0,0); 
thus 0 — —pgx + B, B being an arbitrary constant. 


In passing from the original configuration under zero body force to the reference con- 
figuration under gravity, the free surface moves downwards by the amount 

pgH 2 gH 2 


= r; 0 = 


2K 


Evaluaton of (12) gives 


2c 2 


p° = pg (x - | (H 2 - H t )) +B-B+ = pg(H 2 - x), 


(13) 


(14) 
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which is the correct hydrostatic pressure if 6 « H. If one passes to the incompressible 
limit, c — ♦ oo and 6 — ► 0, but c 2 6 remains fixed and equal to gH 2 / 2. Note that subtracting 
0 eliminates B. The associated displacement field is easily calculated to be 


d° = 



I 


A[ff 2 -(x-/r 2 ) 2 ] o 
0 
0 


+ dr, 


rot > 


(15) 


where is an arbitrary divergence-free rotational motion that satisfies the boundary 
conditions. 


4 Transient Motions 

Subtracting the constitutive relations at the current and reference states we get 

p = — pc 2 V 2 ^> = pc 2 s, (16) 

where s = — V 2 V» is the condensation. Subtracting (12) from (ll) yields 

p - -p(4> - %l>) - ~-[d n } s . (17) 

On equating (16) and (17) we get modified forms of the wave equation that account either 
for nonzero mean boundary surface motions, 

« = ^ (18) 

or nonzero mean dilatation, 

c 2 (V 2 tf>- V 2 ^) = c 2 (V 2 V> + S) =$- 1>. (19) 

The second form follows from —vs = [d n ]s, which is a consequence of the divergence 
theorem. For an incompressible fluid, c — * oo and s = [d n ]s = 0, and either form one 
recovers the Laplace equation V 2 ^> = 0. 

Adjusting the Displacement Potential. 

If the transient displacement potential is modified by a function of time: 

^ = ^ + P(f), (20) 

where ^ is the potential of (6)-(l9), we may chose P(t) so that c 2 \j> = V 2 ^ = — s for any 

A 

t. [P(f) may be found by integrating c 2 iJj — twice in time.] We then recover the 
classical wave equation 

c 2 V 2 ^ = (21) 
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If this adjustment has been made, C(t) vanishes and (17) reduces to 


P = 


-prp. 


( 22 ) 


As an example, consider again the container of Figure 1 in which Hi = 0 for conve- 
nience. At t < 0 the container is in the reference state of rest. At t > 0 it is subjected to 
a prescribed constant-velocity motion of the bottom surface, 

4(*)| =-«*. t> 0, (23) 

positive upwards for a > 0. The unknown free-surface vertical displacement is rj(t), also 
positive upwards. As all quantities become independent of X2 and X3, the governing 
equation is one dimensional: 


*» = £r* + 4 (_ at+ „), 


(24) 


c‘ H 

where primes denote derivatives with respect to Xi = x. The solution for 0 < t < H/c is 


*7 = 0, rl>{x,t) 


aix + ^x 2 + P(t), x < ct, 
P(t) — |ac£ 2 , x > ct. 


(25) 


The adjustment condition gives P(t) = — |acf 2 and consequently rj) = 0 for x > ct. Hence 
for finite c we have = c 2 ^. If c — ♦ 00, the solution approaches the rigid body motion 
77 = at. 
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